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Abstract 


We construct a holographic multi-band superconductor model with each 
complex scalar field in the bulk minimally coupled to a same gauge field. Tak- 
ing into account the back reaction of matter fields on the background geometry 
and focusing on the two band case with two scalar order parameters, we find 
that depending on the strength of the back reaction and the charge ratio of the 
two bulk scalars, five different superconducting phases exist, and three of five 
phases exhibit some region where both orders coexist and are thermodynami- 
cally favored. The other two superconducting phases have only one scalar order. 
The model exhibits rich phase structure and we construct the full diagram for 
the five superconducting phases. Our analysis indicates that the equivalent 
attractive interaction mediated by gravity between the two order parameters 
tends to make the coexistence of two orders much more easy rather than more 
difficult. 
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1 Introduction 


The AdS/CFT correspondence |1, 2, 3] provides us a useful tool to study strongly coupled 
systems holographically in a dual framework through an appropriate gravity theory living 
in a higher dimensional space-time. One of the most studied objects is the so-called 
holographic superconductor. The first holographic model was constructed in refs. [4, 5]. 
Since the condensed field is a scalar field dual to a scalar operator in the field theory side, 
it is a s-wave model. Such holographic approach was also generalized to the p-wave case [6] 
and d-wave case |7, 8]. The basic idea is as follows. Some matter fields, such as gauge 
fields and/or scalar fields, are added into the bulk. One of them considered as “hair” of the 
background solution plays the role of order parameter in the dual boundary system. As 
one tunes some parameters, such as temperature and chemical potential, the background 
solution without hair will become unstable and new stable solution with “hair” appears. 
This process in the bulk mimics the superconducting phase transition in the condensed 
matter theory. 

Such simple holographic setups indeed uncover some basic properties of the real high 
T, superconducting materials. Nevertheless, most of studies existing in the literature have 
been based on a specific setup where the dynamics in the bulk involves only one order 
parameter. It is desirable to generalize the single order parameter case to multi order 
parameter case. 

The high T, superconducting systems, which are thought to be controlled by strongly 
coupled interactions, indeed involve various orders, such as magnetic ordering and su- 
perconductivity, see for example, ref. [9]. The holographic correspondence provides us a 
convenient way to investigate the interaction for these orders by simply introducing dual 
fields in the bulk as well as appropriate couplings among them. The authors of ref. [10] 


studied the competition between different orders in the holographic approach. They de- 
scribed a superconducting order by a charged scalar field and a magnetic order by a neutral 
scalar. It was showed that the appearance of one order inhibits the other one. But when 
the interactions between the bulk fields are repulsive enough both orders can coexist and 
even enhance each other. The authors of ref. [11] focused on the unbalanced holographic 
model by introducing two scalars charged under two U(1) gauge fields respectively. One 
accounts for the electrically superconducting order and the other indicates the electrically 
neutral magnetization. It was shown that the competition and enhancement between the 
two orders correspond to the attractive and repulsive interaction between the two scalar 
fields. In addition, the competition between superconducting and spatially modulated 
phases was exhibited in ref. [12]. 

On the other hand, some new high T, materials such as Magnesium diboride (MgB2) 
and the iron pnictides (LaF'eAsO\_,F,,LiFeAs, Fe ,+.Se, etc.) are characterized by multi 
Fermi surfaces and theoretical investigations have been based on multi-band models ini- 
tiated from refs. [13, 14, 15]. Such multi-band case can be phenomenologically described 
by the Ginzburg-Landau (GL) theory with several complex scalar fields minimally coupled 
to one gauge field. Within this framework, some significant properties of multi-band su- 
perconductors have been revealed, such as fractional flux vortices and vortex bound states 
(see, for example, refs. [16, 17, 18]). One can build a holographic superconductor model 
involving more than one order parameters by straightforwardly generalizing the GL theory 
to the gravity side. Concretely, one can consider a gravity theory with a negative cosmolog- 
ical constant, a gauge field as well as some complex scalar fields. The latter’s condensates 
trigger the superconducting phase transition. In this setup, all scalar fields in the bulk are 
charged under a same U(1) gauge field. Of course, there also exist other ways constructing 
holographic superconductor models with more than one order parameter. Treating three 
bulk scalars as the fundamental representation of SO(3) group, the authors of ref. [19] 
studied a multi-band superconductor model, while regarding two bulk scalar fields as the 
fundamental representation of U(2) group, the authors of ref. [20] investigated a S+/- 
two-band superconductor model, where the AC conductivity exhibits a mid-infrared peak 
which is argued to be related to the inter-band interaction in iron-based superconductors. 

The holographic model involving two competing scalar fields coupled to a single gauge 
field was first discussed in ref. [10], where the bulk scalar fields carry different charges 
and masses, but have no direct interaction between them. The authors of ref. [10] worked 
in the probe limit. Namely, they neglected the back reaction of matter fields on the 
background geometry, the planar Schwarzschild-AdS black hole. It was found that once 
a scalar field condenses, it will hinder the condensation of the other one. ! Depending 
on model parameters, there are two different cases. The first case is that once a scalar 
condenses, the second one will never condense. The second case is more interesting. One 
scalar field first condenses, there the superconducting transition happens. When one lowers 
temperature, the second scalar begins to condense, resulting in a phase where the both 


'The similar phenomenon was also noticed in the case of a non-Abelian model with SU(2) gauge 
symmetry in ref. [21]. 


scalar fields condense. When one further lowers temperature, the first condensation quickly 
goes to zero at a certain temperature and the second condensation goes to a finite value. 
In the superconducting phases, there exist three different cases: first scalar condenses only, 
the second scalar condenses only and both condense. One has to determine which case is 
thermodynamically favored. This can be done by comparing the free energy for each case. 
In addition, the probe limit is expected to be valid only near the critical point. ? In order 
to reach a more persuasive conclusion, one needs to consider the back reaction of matter 
fields on the bulk geometry. It was argued in ref. [10] that once the back reaction of matter 
fields is taken into account, the coexistence region of the two orders will shrink. However, 
our study below shows that the back reaction will reveal more rich phase structure of the 
model. In addition, let us mention here that in a similar model with same mass and charge 
for two scalar fields, it was found [22] that in the probe limit, there is an unexpected gapless 
excitation with quadratic dispersion in the quasinormal mode spectrum. 

We will construct a holographic multi-band superconductor model involving some scalar 
hairs in (3+1) dimensional anti-de Sitter space-time in a rather succinct way by generalizing 
multi-band GL theories to the gravity side. The back reaction of matter fields will be 
considered. As a concrete example, we will consider the two band case. The direct coupling 
between the two scalar fields, in some sense, can be effectively accountable with a simple 
shift of the effective mass of bulk scalars, which would not change the system very much [11]. 
However, the situation is much more complicated by including the back reaction, since the 
effect of gravitational interaction can not be simply considered as only a shift of the effective 
mass for each bulk scalar. Therefore, we will ignore the direct coupling between two scalar 
fields. Of course, one of motivations of this simple setup is to compare with the results 
in [10] and to see clearly the effect of the back reaction of matter fields. Here it is worth 
pointing out that the two scalars have an effective attractive interaction through gravity 
and this attractive interaction becomes stronger and stronger as one increases the strength 
of the back reaction. 

In our study, we can indeed find the coexistence region of two order parameters similar 
to the one studied in the probe limit [10]. Nevertheless, we also find two additional new 
phases in which the two order parameters always coexist once they appear inside the 
superconducting phase. As we increase the strength of the back reaction, our numerical 
calculation reveals that, for suitable parameters, the region of coexistence phases will 
enlarge, rather than shrink, and finally the two order parameters always coexist as one 
increases the back reaction. For both new phases, one of the scalar orders first condenses 
inducing the superconducting phase transition. Then the other order emerges at a lower 
temperature, which triggers a new phase transition within the superconducting phase. The 
similar new phase was also observed in a top down setting in ref. [23]. The calculations 
of the conductivity and free energy uncover that the phase with two order parameters 
coexisting is indeed a superconducting phase and is thermodynamically favored. Adopting 
an eigenvalue method similar to the one in ref. [24], we construct the parameter space for 


?Those studies [10, 11, 20, 19] we mentioned above are carried out in the probe approximation, neglecting 
the effect of matter fields on the background geometry. 


the coexistence phases. The system exhibits rich phase structure. Apart from the normal 
phase, we have totally five different superconducting phases. It is clear that including the 
back reaction is important to complete the full phases of the two band model. The full 
phase diagram for the five superconducting phases is constructed in figure 11. Depending 
on the model parameters, i.e., the strength of the back reaction and the charge ratio 
between the two bulk scalars, each phase is the most thermodynamically stable phase in 
some region of the parameter space. 

The paper is organized as follows. In section 2, we introduce the holographic model, 
obtain the equations of motion of the system and specify the boundary conditions to be 
satisfied. In section 3, we discuss the details of the phase transitions and study thermody- 
namic properties of the system. We calculate the optical conductivity in section 4 to make 
sure that the new phases indeed describe superconducting state. We study, in section 5, the 
parameter space from which one can know the existence range of the coexistence phases, 
and construct the full phase diagram for all superconducting phases. The conclusion and 
further discussions are included in section 6. 


2 The Holographic Model 


In this paper, we will study a holographic superconductor model with N scalar hairs in 
(3 + 1) dimensional anti-de Sitter space-time. The action reads 


= 55 | Pava eee -|V — ier Apl? — mll?) —Vintactl, (1) 
k=1 
where L is the AdS radius, k? = 87G is related to the gravitational constant in the bulk. 
ek and Mg (k = 1,2, ..., N) are the charge and mass of the scalar field eo respectively. The 
term Vintace denotes ilie possible interaction among bulk matter fields. Fi, = V pAr -Vr A, 
is field strength for the U(1) potential A,,. 

We can see that the system depends on the mass mą and charge ex of each scalar field 
Wp. However, one can perform a rescaling of the type A, > = Ais Yk > Eir to set the 
charge of the scalar field Yə to unity. Although we have N scalar fields corresponding to N 
order parameters in the dual boundary theory, the superconducting phase appears as long 
as one of the scalars has a nontrivial configuration. We are interested in the dynamics and 
mutual interaction among different orders. As a concrete example, we limit ourselves to 
the case with N = 2. The model we will study in this paper is described by the following 
action 
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where we have defined D1, = V, -1 A, and D2, = V,,—7tA,. The parameter e> controls 
the strength of the back reaction and e1/e9 is the effective charge of pı or the ratio of two 
scalar charges. 


(2) 


The equations of motion for the action (2) read 
D1“D1 yı — mi = 0, 
D2" D2, — mp2 = 0, 
Vi Fav = i VD Lyi — PDI vi] + ilh Duta — aD ha], 
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(3) 
We would like to find static hairy black hole solutions to these equations, which mimic the 
superconducting phase. 


2.1 The ansatz and equations of motion 


The hairy black hole solution is assumed to take the following metric form 


d? = — f(a? ET + Pde? + dy’), (4) 
f(r) 


together with homogeneous matter fields 


ti=Wlr), p2 = lr), A= ġ(r)dt. (5) 


The horizon r, is determined by f(r) = 0 and the temperature of the black hole is given 
a (ra)/2 

I i 

-r (6) 
One can use the U(1) gauge symmetry to set pı to be real. After using the r component 
of Maxwell’s equations we can also safely choose pə to be real. We will work in the unites 
where L = 1. Then, the independent equations of motion in terms of the above ansatz are 
deduced as follows 
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where a prime denotes the derivative with respect to r. 


6 


2.2 Boundary conditions 


In order to find the solutions for all the five functions F = {q1, Y2, ¢, f, x} one must specify 
suitable boundary conditions at both AdS boundary r — oo and at the horizon r = rp. In 
addition to f(ra) = 0, one must require ¢(r;,) = 0 in order for g” A,A, to be finite at the 
horizon. 

In order to match the asymptotical AdS boundary, the general falloff of the matter and 
metric fields near the boundary r — oo should behave as 


_ Pp o p Piy _ p Wo4 
ọ=upu TARER m= a Aa Ua he aa ee (8) 
€ 
ioe a aa x¥=0+..., 
where Aj, = =?’ ae (i = 1,2). According to the AdS/CFT correspondence, we impose 
w;_ = 0, since we want the U(1) symmetry to be broken spontaneously. Following the 


AdS/CFT dictionary, up to a normalization, the coefficients u, p and Yi, are interpreted 
as chemical potential, charge density and the expectation value of the scalar operator O; 
in the dual field theory, respectively. 

Regularity of the solution at the horizon r = rp requires that all our functions have 
finite values and admit a Taylor series expansion in terms of (r — rp) as 


F= Fira) + F(a) ra) ©) 


By plugging the expansion (9) into (7), one can find that there are five independent parame- 
ters at the horizon {rp, Wi(rn), Yra), O' (Tr), X(ra)}. Note that the equations of motion (7) 
have two scaling symmetries 


eX» eX, ġà, t> AH, (10) 


r= àr, (t,2,y) => A(t, ry) f>Xf, Ad. (11) 


Taking advantage of the two scaling symmetries, we can first set {r, = 1, x(ra) = 0} 
for performing numerics. After solving the coupled differential equations, we should use 
the first symmetry again to satisfy the asymptotic condition y(o0) = 0. Thus we have 
three independent parameters {v1 (rn), Y2(rn), O(n) }, where two of them will be chosen as 
shooting parameters to match the source free conditions, i.e, Y1_ = 0 and %2_ = 0. After 
solving the set of equations, we can obtain the condensate (O;), chemical potential u and 
charge density p by reading off the corresponding coefficients in (8), respectively.’ 

The normal phase in the dual field theory is characterized by the vanishing vacuum 
expectation value of both condensates OQ; and O2, which corresponds to vanishing scalar 


3In our unites, the AdS/CFT dictionary gives (Oj) real = Ot bi and Preal = se In what 
2 2 


follows, we neglect these prefactors and this will not change our conclusions. 


fields yı and Y2 in the bulk. The gravity background describing the normal phase is just 
the Reissner-Nordstr6m-AdS black hole with planar symmetry 


rh i=, 


=al- tr) = alr) sra- pa (02) 


where rp is the black hole horizon and yz is the chemical potential of the black hole. 


3 Thermodynamics and Phase Transition 


As we can see from (2), the two band model is controlled by four model parameters, i.e., 
m?, m3, e2, and e1/e2. We will choose m? = 0 and m} = —2 in this paper. * Following 
the study in [10], one may expect that the model admits three different superconducting 
phases. The first superconducting phase corresponds to the case with Yı # 0 and y2 = 0 
(Phase-I). The second superconducting phase corresponds to the case with wv. 4 0 and 
w, = 0 (Phase-II). The third superconducting phase admits the region where both of the 
scalars condense. 

The first two cases with only one scalar condensing were discussed in ref. [5]. The con- 
densate for a single scalar hair as a function of temperature with a fixed chemical potential 
for such two cases is drawn in figure 1. The condensate curves are very similar to each other 
for different parameters. As we lower temperature, the normal phase becomes unstable to 
developing scalar hair at a certain critical temperature T}. The critical temperature Te will 
decrease if we increase the strength of the back reaction. 
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Figure 1: The condensate as a function of temperature for Phase-I with m? = 0 in the 
left plot and Phase-II with m3 = —2 in the right plot, respectively. The green curves 
correspond to e1/e2 = 2 and e2 = 4, and the purple lines correspond to e1/e2 = 2 and 
€2 = 2: 


4One of the reasons for this choice is for the convenience in numerical calculations. The other is to 
compare our results with those in the probe limit [10]. 
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Figure 2: The condensate as a function of temperature for Phase-C. The red curve is for 
pı, while the blue one is for Ys. We set e;/e2 = 1.95 and eg = 4. 


Except for these two cases with a single scalar condensate, we find other three kinds of 
condensate behaviors where the model admits the coexistence region of two scalar conden- 
sates. The first one is the case uncovered in the probe limit in ref. [10]. This superconduct- 
ing phase denoted by Phase-C is drawn in figure 2. One can see clearly from the figure that 
as we lower temperature, the scalar y1 first condenses at T, where the superconducting 
phase transition happens; when we continue lowering temperature to a certain value, say 
To, the scalar Y% begins to condense, while the condensate of yı decreases, resulting in 
the phase with both orders; if one further lowers temperature, the first condensate quickly 
goes to zero at a temperature, say, 73; when temperature is lower than 73, there exists 
only the condensate of qo. 

The second case is presented in the left plot of figure 3. We can see from the plot 
that the condensate of y first happens at the critical temperature T,; when one lowers 
temperature, the condensate of Y, emerges at a temperature, say To; if one continues 
lowering temperature, the condensate of pə increases while Yı decreases, but the latter 
will never go to zero. In this case, when temperature is less than 7>, both orders always 
coexist. We denote this case by Phase-A. Depending on the back reaction, the inverse is 
also true: the condensate of Y, emerges following the condensate of w2, and then both 
orders are always present. This case is labeled as Phase-B drawn in the right plot of 
figure 3. Whether these three coexistence phases make physical sense depends on whether 
they are thermodynamically favored in their own parameter spaces. Here let us mention 
that in our setup, for a given ex, the critical temperature T, for Phase-B is fixed, while the 
critical temperature T, for Phase-A increases as one increases the value of e1/e2. 
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Figure 3: The condensate as a function of temperature for Phase-A with e;/e2 = 1.95 
and e2 = 2 in the left plot and Phase-B with e1/e2 = 1.9 and e2 = 1.5 in the right plot, 
respectively. The red curves correspond to Yı, and the blue curves correspond to we. 


3.1 Thermodynamics 


Thus we have totally five different superconducting phases in our model. In order to 
determine which phase is the thermodynamically favored phase in some parameter space, 
we should calculate free energy of the system for each phase. We will work in grand 
canonical ensemble in this paper, where the chemical potential is fixed. In gauge/gravity 
duality the grand potential 2 of the boundary thermal state is identified with temperature 
times the on-shell bulk action with Euclidean signature. The Euclidean action must include 
the Gibbons-Hawking boundary term for a well-defined Dirichlet variational principle and 
further a surface counter term for removing divergence 


6 1 4 
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(13) 


where A is the induced metric on the boundary r — oo, and K is the trace of the extrinsic 
curvature. By using of the equations of motion (7) and the asymptotical expansion of 
matter and metric functions near the AdS boundary, the grand potential Q can be expressed 


as 
2K? 


Vo 
with V> = f dxdy. More specifically, for the normal phase shown in (12), one has € = 
=r = me. 

The free energy corresponding to phase-A is drawn in figure 4. Since in this case, 
we has three additional solutions from our equations of motion (7), corresponding to the 
normal phase, Phase-I and Phase-II. In order to make sure whether phase-A is the most 
thermodynamically stable phase, we also plot the free energy of these three phases in the 
same figure. Indeed, the three superconducting phases have lower free energy than the 


= (14) 
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normal phase. Furthermore, Phase-A does have the lowest free energy among the three 
superconducting phases, indicating that once Phase-A appears, it is thermodynamically 
favored. Comparing Phase-B with Phase-A, the only difference is that Y> condenses before 
Wi. Therefore, the free energy of Phase-B also dominates once it appears. 

The free energy corresponding to Phase-C is schematically drawn in figure 5. Similar to 
the above discussion, we also present free energy for other allowable solutions in the same 
figure. As we know, in Phase-C, there is only a narrow window admitting the two orders to 
coexist. Outside this region, it reduces to phases with only a single order. More specifically, 
at the beginning of the superconducting transition, there is no condensate of Y, so Phase- 
C coincides with Phase-I. As Y begins to condense at a lower temperature, the free energy 
of Phase-C becomes lower than Phase-I as well as Phase-II. For a much lower temperature, 
wy, vanishes and Phase-C coincides with Phase-II then. All in all, Phase-C is the most 
thermodynamically favored phase over other possible phases including the normal phase, 
Phase-I and Phase-II. 
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Figure 4: The free energy as a function of temperature for Phase-A, denoted by solid red 
curve. We also plot other possible phases in the same figure. The black dotted line is from 
the normal phase and the green dashed one corresponds to Phase-I. The curve in the insert 
is the difference of free energy between Phase-A and Phase-II. We choose e;/e, = 2 and 
e> = 2. The critical superconducting temperature where Phase-A appears in this case is 
T, ~ 0.0469, which is also the critical temperature of Phase-I. The temperature where Y2 
begins to condense for Phase-A is about 0.0307, and the critical temperature of Phase-II 
is about 0.0451 ju. 
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Figure 5: The typical configuration for the free energy as a function of temperature for 
Phase-C labeled as the solid red curve. In this case, the equations of motion also admit 
other three types of solutions, i.e., the normal phase (dotted black curve), Phase-I (dashed 
green curve) and Phase-II (solid blue curve). There are four special temperature denoted 
as T., Ti, To and T}, corresponding to the critical temperature for the superconducting 
transition, the critical temperature for Phase-II, the temperature at which Y» begins to 
appear in Phase-C, and the temperature where Yı goes to zero in Phase-C, respectively. 
Phase-C and Phase-I coincide as Ty < T < T}, so do Phase-C and Phase-II as T < T}. 


3.2 Superconducting Phase Transition 


As we have seen, for suitable parameters, the solutions with coexisting Yı and Yz appear, 
including Phase-A, Phase-B and Phase-C. Furthermore, such phases are thermodynami- 
cally favored, compared to the superconducting phases with a single order parameter, i.e., 
Phase-I and Phase-II. We now give a further discussion for such three thermodynamically 
stable phases. 

Phase-C is similar to the case found in ref. [10], where the emergency of the scalar y» 
suppresses the condensate of the scalar ~, even although yı condenses before wz. Phase-C 
exhibits three continuous phase transitions. It can be seen clearly from figure 6, where 
the charge density versus temperature is plotted, there exist three particular points at 
which the derivative of the charge density with respect to temperature is discontinuous, 
indicating a second order transition. The one with the highest special temperature is the 
critical point for the superconducting phase transition, while the remaining two points 
inside the superconducting phase. We have checked that the entropy of the black hole 
solution as a function of temperature also exhibits the same discontinuity at these special 
points. 
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Figure 6: The charge density as a function of temperature for Phase-C. The red curve is 

for the normal phase, while the blue one corresponds to the superconducting phase. We 

set e1/e2 = 1.95 and e> = 4. There are three special temperatures at which the derivative 

of charge density with the temperature is discontinuous. 


The condensate behaviors and other thermodynamical properties are similar for Phase- 
A and Phase-B. There exist two transition points in both cases. One is the critical su- 
perconducting phase transition and the other within the superconducting phase, which 
indicates the fact that our system is multi-band. In order to see this clearly, we define the 
total condensate as (O) = (O1,)'/41+ + (O,)1/42+, and draw (O) as a function of temper- 
ature in figure 7. As one lowers temperature, (O) emerges at the critical superconducting 
phase transition point, then at a certain temperature inside the superconducting phase 
it has a sudden increase, where the condensate of the other ~ appears. Such behavior is 
very reminiscent of the one in the real multi-band superconductor, see, for example, the 
figure 3 or figure 12 in ref. [17]. Our discussion can be straightforwardly generalized to 
the case with N > 3 order parameters. For suitable parameters, we believe that the total 
condensate as a function of temperature would exhibit as many as N sudden increases, 
indicating appearance of a new order parameter there. 

It was argued in ref. [10] that as one considers the back reaction, the coexistence region 
for two order parameters in Phase-C will shrink. Surprisedly, comparing the two plots in 
figure 2 and figure 3, one can clearly see that if one increases the strength of the back 
reaction, the region of coexistence phases will enlarge, rather than shrink and the two 
order parameters finally always coexist as the back reaction is strong enough. Phase-A 
and Phase-B are new phases which are not found in the probe limit [10]. In particular, 
the probe limit analysis uncovered that for the chosen masses, i.e., m? = 0 and m3 = —2, 
y, will never condense if Y2 condenses firstly. Nevertheless, our results show that it is not 
the case once the back reaction is taken into account. One can indeed find the Phase-B in 
the right plot of figure 3, where Yı condenses following y2 and then both coexist. 

In some regions of parameter space, the equations of motion do not admit black hole 
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Figure 7: The total condensate as a function of temperature for Phase-A. We set e;/e2 = 2 
and e = 2. The two special points at T, ~ 0.0488 and T œ 0.0298 correspond to the 
superconducting critical point and the point at which Y> begins to emerge in Phase-A, 
respectively. 


solutions corresponding to Phase-A, Phase-B and Phase-C, instead what we can obtain 
is the case with only one scalar condensate, i.e., either Phase-I or Phase-II. For a given 
temperature, the thermodynamically favored phase is obtained by comparing the free en- 
ergy of Phase-I and Phase-I]. One might expect that either Phase-I or Phase-II is always 
thermodynamically favored over the other. However, there might be an interesting possi- 
blity that one of the two phases is thermodynamically favored at the beginning, and as we 
lower temperature, the other phase might become more thermodynamically stable. Thus, 
the physical state is the competition between Phase-I and Phase-II, and the charge density 
as well as the entropy will be discontinuous at a certain point inside the superconducting 
phase, indicating a first order phase transition. In our calculations we have not found such 
a possibility. 


4 Conductivity 


We have found that there exist five superconducting phases in the holographic two-band 
model. In order to ensure the system is indeed in a superconducting state and to see 
whether there are any new phenomena occurring in such a two-band model, we would 
like to calculate the conductivity o. Since we move away from the probe limit, we have 
to consider the fluctuations of A, and gi. Assuming both perturbations have a time 
dependence of the form e~“’, we can obtain the equations of motion for A, and gi, by 
linearizing the full equations of motion (7), which read 
t / 2 2 2 / 9 
Abt (E EA + et EG Ae + Sele ou) =0, (15) 
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fy — “Iie + +o =0. 16 
Its — 79t + ae (16) 


Substituting (16) into (15), we obtain the final equation of motion to calculate the con- 
ductivity 
/ / 2 12 2 e2 
P RE yar y=. 2 E a =: 17 
0+ (GDA, + (G- Sale — Gul + vi) (17) 
Since the conductivity is related to the retarded two-point function of the U(1) current, 
i.e, o = +G*(w, k = 0), we impose the ingoing boundary condition near the horizon 


A, = (r—- rn) ET [ao + ay(r — ra) + aar = ra) +], (18) 
with ag, a1, a> being constants. The gauge field A, near the boundary r — oo falls off as 


AQ) 
Ae = AU ee, (19) 
r 


According to the AdS/CFT dictionary, the retarded Green function can be read as G? = 


(a) l ; ~ 
— from which one can obtain the conductivity 
2 


1 1 A® 
olw) = =G? (w, k = 0) = 


7 WE wA ~ 
The optical conductivity as a function of frequency in the region with two order parameters 
is presented in the left plot of figure 8. For a comparison, we also draw the conductivity for 
the case with only one order in the right plot of figure 8. We can see clearly that the optical 
conductivity in two band model behaves qualitatively similar to the model with only one 
scalar order [5]. In addition, from the Kramers-Kronig relations, one can conclude that 
the real part of the conductivity has a Dirac delta function at w = 0 since the imaginary 
part has a pole, i.e., Im{o(w)] ~ =. 


5 Parameter Space and Phase Diagram 


The calculations of the free energy and conductivity reveal that the phase with coexistence 
of two order parameters is indeed a superconducting phase and is thermodynamically sta- 
ble. More specifically, depending on the model parameters, one gets two kinds of condensate 
behaviors. One is the case that the order parameter emerging at lower temperature sup- 
presses the condensate of the first scalar, the latter will go to zero quickly (Phase-C). The 
other is the case that once the two orders appear, they always coexist, including Phase-A 
and Phase-B. 

It is helpful to construct the parameter space form which one can learn in which region 
the superconducting orders can coexist. However, to numerically find the solution space 
where the two superconducting orders coexist is a little difficult. By our shooting method, 
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Figure 8: The optical conductivity as a function of frequency at temperature T = 0.02734 
for Phase-A (left plot) and Phase-I (right plot), respectively. The red solid line is the real 
part of the conductivity, while the blue dashed line is the imaginary part of the conductivity. 
We choose the parameter e;/e2 = 1.95 and ep = 2 and the critical temperature of Phase-A 
is 7T} œ 0.0469. There is a delta function at the origin for the real part of the conductivity 
in both cases. 


it is more easy to find the single order phase rather than the phase with two orders. 
Nevertheless, one can complete this task by just turning the problem around. Focus on the 
concrete model discussed in our paper, i.e., m? = 0 and m3 = —2, a good starting point is 
to find the critical valve of the ratio e,/e2 such that T is a critical temperature at which 
w, begins to vanish or emerge. At such a temperature, Yı is very small and can be treated 
as a perturbation on the background where only w2 condenses 

' 1 9 wae e? gex 

-v - É -Eriw +n = 9S 


f 2 r 
where {¢ġ, f,x} are functions describing the hairy AdS black hole with only %> non- 
vanishing. 

We demand 7; to be regular at the horizon and to fall off as in (8) near the AdS 
boundary. Then this equation can be considered as an eigenvalue problem with positive 
eigenvalue e7/e3. ° The parameter space is presented in figure 9 and figure 10, where 
different curves correspond to different strengths of back reaction 1/e3. Comparing the 
two figures, we can find that the ratio e,/e2 versus the critical temperature at which the 
solutions with a single scalar %2 become unstable to developing new scalar hair behaves 
differently for large and small back reactions. 

Figure 9 corresponds to the case with large eg, i.e., small back reaction. The critical 
ratio e,/e2 on the curve for a given e> decreases with respect to temperature, and there 
is a one-to-one correspondence between each point on the curve and Phase-C. Concretely, 
in Phase-C, the temperature below which the order y, vanishes is just the temperature 
given here with the corresponding value of e;/e2 that can be read on the curve. One can 


pı = 


pı, (21) 


5A similar method has been adopted in ref. [24]. 
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Figure 9: The parameter space for small back reaction. For each e2, we plot the ratio 
€,/€2 as a function of the critical temperature at which the solutions with a single y 
become unstable to developing new scalar hair. Every point gives the value of e;/e2 and 
the corresponding temperature below which the order yı tends to vanish. The curves from 
up to down correspond to e> = 8.0, 4.0, 2.6, 2.0, respectively. 


further find that each curve has a maximum (e1/e2) maz at T ~ 0 as well as a minimum 
(€1/€2) min at T = Tmar, Which means that Phase-C can appear only when the ratio e;/e2 
is in between them, i.e., e1/e2 € [(€1/€2) mins (€1/€2) max: For the case e1/e2 > (€1/€2) mars 
Phase-A appears and dominates the system. If we decrease the value e;/e2, Y2 will condense 
before yı. Our eigenvalue method tells us that the superconducting background with w» 
hair only does not become unstable under the perturbation imposed by %1, which means 
yy, can never appear inside the condensate of w2. This is just the case denoted as Phase-II. 
Nevertheless, the case with Yı 4 0 and Y3 = 0 labeled as Phase-I is also the possible phase. 
Thus, we have two simple phases with only a single order, i.e., Phase-I and Phase-II when 
e1/e2 < (€1/€2),,,;,. For a given temperature, the thermodynamically favored state is the 
consequence of the competition between Phase-I and Phase-II. As we discussed earlier, one 
of the two phases would always be thermodynamically preferred to the other. And the 
other possibility is that one of the two phases is firstly thermodynamically favored, then 
the other phase becomes more thermodynamically favored, which indicates a first order 
phase transition. However, our numerical calculation rules out the later possibility and 
the thermodynamically favored state is Phase-II. As one can see, as the strength of the 
back reaction increases, the value of (e€;/e2),,,;, and (€1/€2),,4, for each curve increases and 
decreases, respectively, which indicates that the attractive interaction mediated by gravity 
reduces the parameter space of Phase-C, but enlargers the parameter space of Phase-A. 
Let us now move to the case with large back reaction, which is presented in figure 10. In 
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Figure 10: The parameter space for large back reaction. For each e2, we plot the ratio e1/e2 
with respect to the critical temperature at which the solutions with a single Y> become 
unstable to developing hair yı. Each point gives the value of e;/e2 and the corresponding 
temperature below which the order y, begins to appear. From the up to down in the 
leftmost the curves correspond to e = 1.5, 1.3, 1.1, respectively. 


contrary to the small back reaction case, for a given e2, the critical ratio e;/e2 on the curve 
increases with respect to temperature and has a maximum (Tmar, (€1/€2) maz) and a mini- 
mum (T © 0, (€1/€2) nin). In figure 9, the point located in the range [(€1/€2) min, (€1/€2) maz] 
is related to Phase-C, while the point in figure 10 here corresponds to Phase-B. For 
e1/e2 < (€1/€2)min and e1/e2 > (€1/€2) max, One can obtain phases with only one order 
parameter. The physical state might be in Phase-I, Phase-II or the combination of both. 
Nevertheless, by comparing the free energy between Phase-I and Phase-II, we find that 
Phase-I is thermodynamically favored for eı/e> > (€1/€2) maz and Phase-II is thermody- 
namically favored for e;/e2 < (€1/€2),,;,. We does not find any suitable values for e;/e2 
and e2 for which one can have a phase similar to Phase-C. In contrary to the small back 
reaction case, we find that, as one strengthens the back reaction, (e1/e2),,,;, decreases and 
(€1/€2) ma, Increases, respectively. Therefore, the attractive interaction between the two 
bulk scalar fields mediated by gravity enhances the parameter space of Phase-B. Compar- 


ing the two figures, it is clear that there must be a certain value of e> denoted by es" 


at which the curve becomes parallel to horizontal axis. At e§”"°"’, (e;/e2),,,, has the same 
value of (€;/€2),,,;,, thus there is no region for the Phase-C in figure 9 as well as Phase-B 
in figure 10 to survive. In our model discussed in this paper e§"“°" ~ 1.762. 

To summarize, apart from the normal phase, we obtain as many as five superconducting 
phases in our model. Depending on the model parameters e;/e2 and e, each phase can be 


most thermodynamically stable in some region of parameter space. The full phase diagram 
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Figure 11: The full phase diagram for the five superconducting phases. Depending on e;/e2 
and es, the phase diagram is divided into five parts. We label the most thermodynamically 
favored phase in each part. 


for the five superconducting phases is constructed in figure 11 which is divided into five 
parts by five boundaries. The corresponding phase we named in each region is the most 
thermodynamically stable phase. The solid red line between Phase-C and Phase-II and the 
dashed red line between Phase-B and Phase-I is obtained by finding the value of e;/e2 at 
which the two critical temperatures for the phases with only %ı or Ys condensate are the 
same for each given e2. The solid green curve separating Phase-C from Phase-A and the 
dashed green curve separating Phase-B from Phase-II correspond to the leftmost points in 
figure 9 and figure 10. The vertical blue line dividing Phase-A from Phase-I is given by 
the critical line at egea. We see from figure 11 that when the value of ez increases, the 
upper bound of Phase-C (solid green curve) increases, but finally it arrives at a finite one 
around 2.164. In contrast, the lower bound of Phase-C (solid red curve) decreases, and 
finally it arrives at a finite one around 1.867. ê 

As we increase the strength of the back reaction, the region for Phase-C with the co- 
existing behavior of two order parameters only in a narrow window is gradually forced to 
shrink and finally vanishes at e$"“*°“, while the regions for Phase-A and Phase-B where 
both order parameters always present enlarge. In this sense, we can conclude that the 
gravity which provides an equivalent attractive interaction between the holographic order 
parameters tends to make the coexistence of two orders much more easy rather than more 
difficult. The other way to see this is to compare different curves for given ratio e;/e2. We 


Due to the lake of the numerical control for very large back reaction, we can not fix the phase boundaries 


for very small e2. However, our current results indicate that the boundary of Phase-B tends to expand as 
one increases the strength of the back reaction. 
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can find from figure 9 that increasing the strength of the back reaction lowers the temper- 
ature at which Yı vanishes. Thus the coexistence region in Phase-C can survive at much 
lower temperature. On the other hand, in the case with not very small e;/e in figure 10, 
we can see that increasing the strength of the back reaction uplifts the temperature where 
ə background becomes unstable to forming yı hair, resulting in a much wider region 
where both order parameters can coexist. 


6 Conclusion and further discussions 


We have investigated a holographic superconductor model with more than one order pa- 
rameter in four dimensions, where each complex scalar field in the bulk is minimally coupled 
to a same U(1) gauge field. This can be interpreted as a holographic multi-band supercon- 
ductor model. We have studied the interplay among different macroscopic orders in this 
model, where the back reaction of matters on the background geometry has been taken 
into account. 

Concretely, we have discussed the two-band case with mass squares m? = 0 and m3 = 
—2 for two bulk scalar fields yw, and we, respectively. Depending on the strength of the 
back reaction 1/e3 and the relative charge ratio e,/e2 of the two scalar fields, the model 
admits as many as five different superconducting phases. Three of them, denoted by Phase- 
A, Phase-B and Phase-C, exhibit the coexistence region of two order parameters. More 
specifically, for Phase-C, as we lower temperature, the second scalar Y> condenses following 
pı will completely suppress the condensate of the first order, i.e., Yı will go to zero finally. 
The condensate behaviors in Phase-A and Phase-B are similar. One of the two orders 
condenses first, and once the other begins to condense, both always coexist. Other two 
superconducting phases, labeled as Phase-I and Phase-II, have only one order parameter. 
For details, see figure 1, figure 2 and figure 3, respectively. 

For given parameters, i.e., e,/e2 and e2, Phase-A, Phase-B and Phase-C can not appear 
at the same time, while Phase-I and Phase-II are always allowable. Therefore, once we 
have a phase with coexisting orders, we can also have the two phases which have only one 
order. We have calculated the free energy for each possible phase and found that the three 
phases with coexisting orders are thermodynamically favored than Phase-I and Phase-II. 
For each ez, we plot the ratio e,/e2 as a function of the temperature at which the solutions 
with only p condensate become unstable to developing new scalar hair Y, in figure 9 
and figure 10. The behavior of e;/e2 versus temperature is exactly opposite for small and 
large back reaction cases. The model shows rich phases structure and we construct the 
full phase diagram in figure 11. From such three figures, we can conclude that the gravity 
providing an equivalent attractive interaction between the holographic order parameters 
in some sense will enhance the coexistence region and tend to make the coexistence of two 
orders much more easy rather than more difficult. It also implies that to reveal a full phase 
structure for a holographic model, the back reaction of matter fields in the bulk has to be 
taken into account. 

In this paper, we have considered in some detail the case with m? = 0 and m3 = —2. We 
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expect that the whole picture will not be qualitatively changed for other choices of mass 
parameters if one always regards the scalar with larger mass as w, and the scalar with 
smaller mass as W2. A special case is that the two scalars have equal mass. In this case, 
both condensates can appear at the same temperature for some parameters, for example, 
e1/e2 = 1 even with no direct coupling between two scalar fields. One may worry that in the 
region which only admits the solutions corresponding to Phase-I or Phase-II, there might 
be a first order phase transition between Phase-I and Phase-II by comparing their free 
energy. We have checked a lot of parameter values and found that either Phase-I or Phase- 
II is always dominant. Nevertheless, even although such a first order transition might exist, 
the full phase diagram does not change. It is interesting to understand whether this first 
order phase transition could appear or not for other possible values of mass. In addition, 
in our discussions, we have turned off the direct coupling between two bulk scalars, it is 
also interesting to study the case with a direct interaction between the two scalars. Of 
course, our study can also be straightforwardly generalized to the case with more than two 
order parameters. In that case we expect much richer phase structure will appear. 

According to the symmetry of the macroscopic wave function or condensate of Cooper 
pairs in the real superconducting materials, the superconductor can be classified by s-wave, 
p-wave, d-wave and so on. The holographic s-wave, p-wave and d-wave superconductor 
models are already available in the literature, it is therefore quite interesting to study 
the holographic models with superconducting order parameters with different spins. Fur- 
thermore, we may study the holographic multi-band superconductor model by including 
the effect of lattice [24]. In our present study, the gravity background is chosen to be 
black holes, which corresponds to the superconductor/conductor phase transition at finite 
temperature. Note that the holographic superconductor/insulator phase transition at zero 
temperature has been studied in ref. [25] by choosing the so called AdS soliton background 
with one spatial direction compactified to a circle and the complete phase diagram for such 
a system was constructed in ref. [26]. Certainly it is very interesting to investigate the holo- 
graphic superconductor/insulator phase transition with more than one order parameter. 
We leave these issues for further study. 
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